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1 . Introduction 

This paper describes a mathematical analysis for 
determining the value of the substitution error of a 
bolometer with a Wollaston-wire element (barretter). 
The analysis reflects all significant nonlinearities in 
the heat flow, including some not covered before, and 
includes all appreciable heat transport mechanisms 
simultaneously. The microwave current distribution 
has been measured in a scaled model and found to be 
approximately sinusoidal in shape with a current 
maximum in the middle of the wire. Input data are 
chosen so as to represent the operating bolometer 
characteristics as nearly as possible. 

The values of substitution error thus obtained, in 
conjunction with efficiency data obtained by micro- 
wave techniques [l], 1 will be very useful in extending 
power meter calibrations to frequency ranges where 
extremely accurate microcalorimeters are not avail- 
able. A method of measuring substitution error in- 
volves two distinct measurement techniques [1, 2], 
and since agreement is quite good between the value 
of substitution error thus measured and the value 
calculated from this analysis, the two provide a cross- 
check for each other. 



2. Background 

In the measurement of microwave and millimeter- 
wave power (hereafter referred to as rf power) by 
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substitution techniques, any substitution error is of 
great concern. A common measurement technique is 
to replace a known amount of d-c power with an 
unknown amount of rf power in a bolometer. Un- 
fortunately, the different current distributions generate 
different temperature fields which give the bolometer 
element slightly different values of total resistance 
for equal amounts of power. For balanced-bridge 
methods the resistance is maintained constant, which 
causes some non-equivalence of power. In this case 
Carlin and Sucher [3] defined substitution error as 

where W xi is the rf power dissipated in the barretter 
element and JF S ub is the substituted d-c power which 
yields the same value of barretter resistance. Although 
a bolometer can be used to measure relative powers 
without calibration, it must be calibrated to measure 
absolute power. An accurate method of calibration 
is by the use of a microcalorimeter [2], in which the 
ratio of J^sub to the net rf power flowing across an 
arbitrary plane J^net into the bolometer is measured. 
This quantity is called the effective efficiency, 
t)e— ^sub/^net- Another calibration method (for 
barretters only) is the impedance technique [1], in 
which the ratio of W xi to J^net is measured. This quan- 
tity is called the efficiency, r\= JV rf IJV NKT . By replacing 
J^sub and W x t in the equation for E, an equivalent form 
is given by 



E = 
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This provides a way of measuring substitution error, 
but with the substantial limitation that in typical cases 
this necessitates subtracting two nearly equal numbers, 
and the relative uncertainty of the result is very large. 
Some results, comparisons, and comments will be 
given later. 

Another equation derived from the previous ones 
is given by 



"NET 



r sUB 

T?(l+£) 



This form provides a method of measuring the net 
power from W SVB , V, an d E. The substituted d-c power 
is measured relatively easily, efficiency can be meas- 
ured by the impedance technique (as mentioned 
previously), and now with the results of the analysis 
of this paper, substitution error can be calculated. 
For frequency ranges where microcalorimeters are 
not available, or for organizations to which they are 
not accessible, this provides an accurate means of 
calibrating barretters. 

Previous work on determination of substitution 
error has been done by Carlin and Sucher [3], Weber 
[4], and Bleaney [5]. However, none of the previous 
workers attempted to solve the total heat-flow problem 
or to take into account all of the nonlinearities. Carlin 
and Sucher have determined an upper bound of sub- 
stitution error from a combination of solutions to 
parts of the total problem. Weber has handled one 
nonlinearity in a solution to part of the problem. 



3. Description of the Barretter 

The actual barretter as shown in figure 1 can be 
represented by the model shown in figure 3. Partial 
justification for this is based on results of a model 
that used electric potential to represent temperature 
potential (fig. 2). These results showed that flush sup- 
ports in a model yield negligible difference in potential 
from other common support structures. This is im- 
portant, not only because flush supports give the 
simplest mathematical model, but also because there 
are large variations in the shapes and dimensions of 
support structures, even for the same models by the 
same manufacturer (fig. 4). 

The rest of the justification for using the model 
shown in figure 3 becomes clear in the discussion of 
operation and different heat transport mechanisms 
that follows. 

The barretter is a small-diameter, platinum wire of 
radius a and length / (fig. 3). It is supported at each 
end by a silver wire of larger diameter, and is normally 
immersed in air. Possible external, forced-convection 
currents are eliminated by a plastic shield of radius b 
which is nominally transparent to the electromagnetic 
fields. 

For power measurements, this fine wire is usually 
heated with a d-c current to a particular value of 
resistance in the absence of rf fields, and then the d-c 
current is reduced to maintain the previous value of 
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(Fine Nichrome Wire 

Electrolyte 

Metal Rod, Diameter 
to Scale of 
Silver Support 



Metal Tub 



-Metal End Plate 
FIGURE 2. Bath for temperature-potential analogy. 

= 0° 




Ambient Temperature is 0° 
in this Case 

- k 20> ^21 




FIGURE 3. Model used to represent barretter. 



Hrrr 



fm 



FIGURE 4. Various common support structures. 
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resistance when the rf field of frequency / is applied. 
The substituted d-c power is (/f — !:,)/{, where h is the 
d-c current without rf power, / 2 is the d-c current with 
rf power, and R is the value of resistance to which the 
barretter is biased. The substitution error is caused 
by the difference in temperature fields between the 
two heating conditions. The solution of these two tem- 
perature fields will then require the additional data 
concerning thermal properties of the wire and the 
surrounding fluid, usually air. The thermal conduc- 
tivity of the wire, with its temperature coefficient, 
may be represented by ki = k\o-\- k\\T. Similarly the 
thermal conductivity of the fluid may be represented 
by ko = A~2o + A 2 1 7 1 . The resistivity of the wire and the 
corresponding temperature coefficient may be repre- 
sented by p = p 4- p i T. 

The previously mentioned parameters could all 
be input data; however, as is mentioned later, it is 
more judicious and accurate to arrive at some of these 
parameters from other measurements taken under 
operating conditions. 

All heat-transfer mechanisms that can possibly act 
in a barretter will have to be considered and evaluated. 
Of conduction, convection and radiation, only conduc- 
tion turns out to be significant. Radiation is negligible 
since surface area of the wire is so very small. Power 
radiated is given by 

p=aSF(T$-Tl), 

where p is power in watts, 

a is the Stefan-Boltzman constant, 
S is the surface area of the fine wire, 
F is a surface emissivity factor and the 7"s are 
absolute temperatures in °K. 

The diameter of the wire is of the order of magnitude 
10 -4 cm, and typical calculations show the radiated 
power to be less than 0.01 percent of the total power, 
even if the wire were at its melting temperature. 
Convection may also be shown to be negligible (0.01%) 
by evaluation of the Grashof number which is some- 
what proportional to convective heat transport. This 
may also be confirmed by varying the orientation while 
measuring d-c power. The Grashof number is 
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where /3 is the cubical expansion coefficient, 
g is the gravitational constant, 
p is the density, 
L is the width of the gas space inside the plastic 

convection shield, 
17 is viscosity, and 
T's are the absolute temperatures. 

The Grashof number thus calculated for typical 
bolometers is less than 10; if it is less than 1000, 
convection may be considered negligible. 

The only remaining heat transport mechanism, 
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FIGURE 5. 8 X 16 in Scale model of waveguide bolom- 
eter for measuring rf current distribution in 
element. 

conduction, can be through both the fine wire and the 
air (or other fluid). The thermal conductivities of both 
media are functions of temperature over the rather 
wide temperature range for normal operation. 

It is essential to consider the current distribution 
in the wire. The d-c current distribution is considered 
uniform both axially and radially. The rf current is 
essentially uniform radially, even to frequencies as 
high as 100 GHz (< 10% variation), primarily due to the 
smallness of the wire. The axial rf current distribution 
is of greater uncertainty, and has been investigated 
extensively with a scale model (fig. 5). Details are 
discussed in the following three paragraphs. 

The measured values indicate a current maximum 
near the middle of the wire and in general a sinusoidal 
or parabolic shape. There is a frequency dependence 
on the amount of curvature. There is great variation 
in current magnitude but only slight variation in shape 
with respect to relative position of the shorting end- 
plate and with respect to position of tuning screws. 

Some variation in the location of the current maxi- 
mum may be introduced if the shorting end plate is 
tilted only a few degrees from perpendicular. The 
shape and size of the wire support structure cause 
variations in the current shape, but to a large extent, 
only insofar as the wire length itself is varied. An 
extreme case not encountered in most bolometers is 
when there is no support structure — the wire supported 
only by the wide walls. In this case the current dis- 
tribution is nominally uniform, with great susceptibility 
to asymmetry introduced by varying angles in the 
shorting end plate. 

The measured values seem to follow a sinusoidal 
distribution with wavelength greater than free-space 
wavelength and with some asymmetry (< 30° for com- 
mon support structures). The sinusoidal distribution 
has been used by previous workers [3, 4, 5], but effects 
of asymmetry have not been previously considered. 
Further investigation is being considered, but for 
present calculations rf current is assumed to be sinus- 
oidal with a free-space wavelength. 
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There are several pitfalls in determining three of 
the input parameters a, p , and pi. 

First, the radius a or diameter d cannot be 
measured directly with a microscope without 
bringing the objective lens destructively close to 
the barretter if sufficient magnification is used. 
Second, the small diameter of the wire causes 
po to vary from the bulk properties of the metal 
[6] since the wire diameter approaches the mean 
free path of an electron; pi does not vary this 
way [7, 8]. 

Third, a set of three equations with a, p , pi 

unknown is homogeneous and the determinant 

vanishes; it is possible to solve only for two of the 

unknowns in terms of the third. 

The second circumstance indicates that if a value 

for one of these three unknowns must be assumed, 

assuming a value for pi causes the least uncertainty. 

The equation 

Ri7Ta 2 — pol = piTil, £ = 1,2 

can then be solved for a and po if two sets of R and T 
are given; / can be measured accurately. One of the 
R and T sets is for room temperature while the other 
is for operating temperature. The average operating 
temperature is calculated by the computer program, 
as mentioned later. 

4. Analysis 

In this section, we give a description of the analysis 
used to obtain the substitution error and of the approxi- 
mations involved. Details are given in appendix A, to 
which reference is frequently made. 

Referring to the model shown in figure 3, the basic 
mathematical problem is to solve two distinct cases 
of the heat flow equation 



V-(A* + Ai7)Vr=S(r,*, T) 



(1) 



for the temperature field 7(r, x), where (/c , k x ) are 
conductivity coefficients (different for r<a, r> a) 
and S is the local power dissipation. (A— 2, A — 3) 

In the first case, the temperature distribution Tj is 
generated by the source S x (A-l) due to dissipation of 
the d-c current J>\ flowing uniformly in the wire with 
resistivity linearly dependent on T. 

In the second case, the temperature distribution 
T u is generated by the source S u (A— 4) due to dissi- 
pation of a d-c current J^ 2 , as above, and also due to 
the time-average dissipation of an rf current (A-9a): 



J>r{%, t) = l #Roz{x) COS (lit, 



(2) 



where z{x), the current profile, is given, and the 
amplitude J'no is an unknown constant. 

Boundary conditions on the temperature field T(i\ x) 
in each case are: 



T, r (0,x)=T(b,x) = 0, 



(0 



0; 



T(r,0)=T(r,t) = 0, 

{T(r, x), T, r (r, x)} 



(0 = 



b); 



(3) 



continuous across r = a, (0 ^ x ^ /). 

The rf current amplitude J> RQ is to be determined 
by requiring that, given J> x and J^ 2 , the total wire re- 
sistance R hot be the same in each case (A-20): 



R 



hot" 



1 f l 

— - dx(p -f pi7\0, x)). 

Td~ Jo 



(4) 



Because of the nonlinearities, the problem posed 
above is not tractable, and several simplifying assump- 
tions are introduced which lead to a linear system of 
equations involving a single parameter, the mean 
axial temperature: 



h/: 



dxT^O, x). 



(5) 



The linear system is solved iteratively and converges 
rapidly. 

Before making the simplifying assumptions, it is 
useful to introduce the "quasi-temperatures" (A-5): 



G(r, x) = T Y {r, x) 



//(r, x) = T n (r 9 x) 



1 k\ 
1 + 2 T T n( r ' x ) 



g{r, x) = H(r, x) — G(r, x). 



(6) 



These quantities simplify the heat flow equations, 
since 






(7) 



The interface conditions at r=a are, of course, more 
complicated: 



IcqG, r , G ( 1 + t 



h 



(8) 



are continuous across r=a, (0 * 
conditions on H. Here r(x) is 
(A-7b): 



s x ^ /) , with identical 
the inverse function 



G ]-H 



H)- 



(9) 



The three basic simplifying assumptions that are 
made are the following: 

First, the difference function g is assumed small: 



\g(r 9 x)\<\G(r,x)\. 



(10) 



This approximation appears to be adequately justified 
by the fact that computed values indicate \g\ < 0.15|G| 
at worst, while the axial average g=0. Neglecting 
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^(|g| 2 ), we obtain equations of the form (A-8a, A-ll): 
V 2 C-r-^iG[l + ^i(G)] = ^ 2 , 

V>g + A :i g[l+<p->(G)] = e(x) 



[po + piG(l-h 93(C))], 



(11) 



where the A\ are constants (different for r > a, r < a), 
the functions ^(G) depend on r [77 Gj, and for cases 

studied, we have again \<Pi(G) \ < 0.15. 
The function e(x) depends on the difference in current 
flows in the d-c and rf cases. The second (nonlinear) 
interface condition is (A- 14) 



|G(1 + <p 4 (G))\ a a -=0, 



(12) 



\g(l + (ps(G 



= 0, 



with <pi(G) as above. 

The second simplifying assumption is that for 
^ r < a, we may put 



<Pi(G) = <pi(G(0, x)) (13) 

using the axial average, while for r>a, we may put 



(p&G) = <ptG*) 



(14) 



where G* satisfies the interface condition in the mean 
G*(1 + ^4(C*))=G(1 + 94(G)). (15) 

This approximation makes use of the expected flatness 
of the temperature distribution in x 9 and is clearly 
better for large l/b. Numerical tests indicated the 
results were very insensitive to the choice of G; a 5 per- 
cent change in (7 gave a 0.06 percent change in E in a 
typical case. 

The third simplifying assumption, used indirectly 
above, is that we may put 



G(r, x) = G(0, x) 



(16) 



in the source term for g, (r ^ a). This is readily justi- 
fied by the extremely small radial temperature drop 
in the wire. 

With these simplifications, the linear equations with 
parameter (J in (g, G) (A-12, A-13, A-17, A-18) are 
readily solved in terms of cylinder functions. With 
the same simplifying assumptions, we find the balanced 
bridge requirement (A-21): 



Sl0,x) = 



(17) 



which determines the unknown rf current amplitude 
JPro implicit in e(x). — 

After iteration to determine the parameter G, the 
required power levels f^suB ana> ^rf> tne substitution 
error E (A-27), and the temperature distributions 
r,(0, x) and A7 1 = r„(0, x) - r,(0, x) (A-34a, A-34b) are 



readily found by numerical integrations involving 
G(0, x) and the rf current profile z(x). 

Details of the computer programs are available from 
the authors. Further details for the sinusoidal rf 
current distribution are given in (A— 37, A-38, A-39). 

5. Uncertainties 

With given input data, the analysis and computation 
described in the previous section is believed to yield 
a substitution error which is correct for the cases 
studied to within a few percent, though it has not 
been possible to construct an error bound. 

An evaluation of the total differential reflecting the 
uncertainties of all the input data would seem to offer 
a limit of uncertainty to the calculated value of substi- 
tution error, so that one concerned with applying this 
technique would have sufficient justification for doing 
so. This figure is nominally 20 to 25 percent, which is 
quite adequate for the typically small values of sub- 
stitution error. 

Unfortunately the shape of the rf current distribu- 
tion (assumed sinusoidal with a free space wavelength 
and with a current maximum in the middle) is known 
no better than it can be measured with a scale model. 
Uncertainties due to variation of wavelength and shift 
in position of current maximum are partially reflected 
in this total differential; phase shift variations cause 
only 0.2-0.3 percent decrease in substitution error, 
and are included in the above uncertainty, but the 
uncertainty as to wavelength could cause 10 to 40 
percent change either way and is in addition to the 
above uncertainty. Even this is just a best-judgment 
estimate. 

This degradation could be substantially reduced if 
a sound theoretical analysis of the current distribution 
in the barretter were available. Since none is, further 
investigation is certainly in order. 

As to the input data, R cold and R hoi may be measured 
quite accurately, while only two of the four quantities 
p () , a, R co \a and R hoi are independent. First, the follow- 
ing relation must hold: 

R, .=^4- (18) 

A second relation results as follows: Given the param- 
eters a and Xi={R coUU pu h h b, /, ^f, f\, /, tan <p} 
the temperature field T(x, o) is determined in the wire, 
as is the heated wire resistance R 



R 



I 



(po + piT(x,0)) 



iTcr 
- R(a: xt) 



and a substitution error is determined 
E = g(a;xd. 

An iteration is made to find a* such that 
R(a*;xi) = R hot 



(19) 



(20) 



(21) 
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which is the required second relation. Then also 



(22) 



is the proper substitution error. 

These forms are also used to determine the effect 
on E of small uncertainties in the data Xu From small 
changes in (%u a), we compute the partial derivative 
approximations: 



a^ dR dW dR 

da' da' dxf d Xi 



With the bridge balanced 
dR 



v^hile 



dR 
8« hot =— 8a+^8** = 



ot =—- oa + — oxi. 

da dxt 



(23) 



(24) 



(25) 



Then 



8E = 8xi 



dxi 



dR 

d^dXi 

da dR 

da 



(26) 



gives the error in E due to small input error dxi. 

6. Results, Comparisons, and Evaluation 

Values of substitution error have been calculated 
for several common barretters and where possible 
comparisons are made with some measured values 
(see table 1). Although the agreement is reassuring, 
the uncertainty in the measurement makes the com- 
parison inconclusive in itself. 

A number of variations of input data lead to several 
interesting observations concerning the calculated 
values of substitution error. 

Based on calculated results which use sinusoidal 
current distribution with free-space wavelength, 
substitution error varies approximately as the square 
of the frequency. 

The effect of asymmetry in the current distribution 
(induced by asymmetries in the physical structure) is 
to slightly reduce the substitution error, typically of 
the order of magnitude of 0.1 to 0.3 percent. 

The use of helium, a gas with higher thermal con- 
ductivity than air, greatly increases the power han- 
dling capacity of the barretter. It also increases the 
heat "coupling" to the surroundings and thus reduces 
the substitution error approximately 10 percent. The 
same effect is obtained (but without increasing the 
power handling capacity) by increasing the bias cur- 



rent (e.g., so that the resistance increases from 200 
to 267 ft). 

A change in rf power level causes virtually no change 
in substitution error. 

The nonlinearity due to the temperature dependence 
of the thermal conductivity of the air causes an ap- 
proximate 15 percent increase in substitution error 
over what would have occurred if that nonlinearity 
had not been present. The temperature dependence 
of the thermal conductivity of the metallic element 
causes only slight variation (~ 1%) in substitution 
error as only about 5 percent of the heat escapes the 
element by metallic conduction. 



Table 1. 



Comparison of measured and calculated values of 
substitution error 



Freq 
(GHz) 


NBS 

barretter 
unit 


Effective 

efficiency 

iVe) 


Efficiency 

(v) 


No. of 
meas. 


Measured 
substitu- 
tion error 

£ =^-? 


Calculated 
substitu- 
tion error 


8.2 
12.2 


18 
23 


0.9920 ±0.002 
.9709 ±0.002 


0.9906 ±0.005 
.9644 ±0.005 


12 
3 


0.0014 
.0067 


0.0020 
.0056 



It is also of interest to compare results with those 
of previous workers, in particular Carlin and Sucher 
[3]. Although they obtained an upper bound for sub- 
stitution error, and this upper bound is typically 
20—25 percent higher than our calculated results, a 
closer examination of their work points to a flaw in 
one of their conclusions, namely that nonlinearities 
have little effect in these cases. The margin of safety 
in their upper limit comes from the fact that they solve 
two related linear problems, calculate the substitution 
error for each case and add the results together to 
obtain a safe upper limit. It happens that both of these 
related problems resemble the real problem, their 
equation (15), quite closely, and the results of their 
equation (15) should be a good solution to the linear 
approximation to the real problem. The corresponding 
results of equation (15) are about 40 percent lower than 
our results in most typical cases; this indicates the 
magnitude of the error made in assuming that non- 
linearities can be neglected. 



Appendix A. Analysis 

A wire of radius a and length / terminates on me- 
tallic surfaces at z = and z = /, the whole enclosed 
in a shield of inner radius b, as shown in figure 3. The 
wire is heated by current sources uniform in 
r, (0 ^ r^ a), and the temperature T above ambient 
is assumed to vanish on the shield r=6, (O^x^l) 
and on the end surfaces z = (0, /) , (0 ^ r ^ b) . 

Let Ti be the temperature field in the d-c-only case, 
in which the source strength Si in the wire due to d-c 
current $\ is 



s,=>, 2 (po+p 1 r 1 ) 



(A-l) 
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where (p , pi) = (po/A 2 , pi/A 2 ) and A is the cross- 
sectional area of the wire. 
The heat flow vector is 

ft = -(*n + far I )V7 , I (A-2) 

and governing equation for heat flow is 

V-ft=S, (A-3) 

with subscript i referring to ^ r < a and a < r ^ b 

(fig. i). 

Let jf n be the temperature field in the rf-substitution 
case. The source strength S u (x) is taken to be a time 
average due to a d-c current J^ 2 and an rf current 
^r(x, t): 

Sn = ([S 2 + S R (x, t)] 2 )(po + piT n ). (A-4) 

Equations analogous to (A-2, A-3) hold for the heat 
flow in this case. 

The fields (7*,, T u ) and the heat flow vectors (q u q u ) 
are to be continuous across r= a(0 ^ x ^ /) . 

Let us introduce new variables 



l+|-n7i 



H(r, x) = T u ( 1 + 2 nTu) (A-5) 

g(r,x)mH{r,x)-G(r,x), 

/here 



H — /'/i //«'(" 



(A-6) 



We shall suppose the form of the rf-wave is given, 
but its amplitude J> m is determined (for given J>\, J> 2 ) 
to balance the bridge after substitution 

J R (x, t) =Sroz(x) cos co/. (A-9a) 

Then taking time average 

(lf% + Sm*(x) cos co/] 2 )-^ + y.%z i (*) 



Let ( 



so that 

1 



denote the average on (0, /) and define 



zi(x) 



z zHx) 
z 2 (x) 



-1, Zl (x)=0, (A-9b) 



^:^2_lI 



M^^%^{x)^l^^ m zHx){\ + z x {x)). (A-9c) 

In the long wave-length limit, z 2 {x) is a constant, 

1 



so that the definitions differ in the regions r ^ a. The 
relations in (A-5) may be inverted to give, e.g.: 



7-,=G(l + T(r<G)), 



(A-7a) 



r(x) 



1 
~2* 



zi(x) = 0, so that if ^l+j^ z 2 U)=^?, the sources 

S| and S„ are identical, so are the fields 7',, 7",,, and 
there is no substitution error. We shall assume then 
that 

l z Jh^)M-i-J?l + I\ (A-9d) 

and thus that 

S II =Ul+U?-^ + / 2 )(l+2,U))](p ( , + p 1 7 , „). 

(A-9e) 



where PI{J\ — J\) and z x (x) are assumed small com- 
pared with unity. We must expect g(r, x) to he 
small also whenI 2 /(./f-J^|) andz,(x) are small, and 
we linearize (A-8b) accordingly in g: 

k ia &g = - 8« {.?\Bg +BCg+AC], 



l+VH* 



In terms of the Laplacian 



(A-7b) where 



ld_ d_ d\ 

r dr dr dx 2 



~ »- x r r z r ' ;w2' 



/4 = po+/5iC(l-fr(r 1 G)) 

fi = Pi(l + T(r 1 G+r 1 Gr , (r 1 G)) 
C=/2+(/f_/|) Zl ( A; ) + ^(/2 2l ( JC )) > (A _ 10) 



kioJz 



Sl6n[fio + piG(l + T(r t G))] 

8i,<(/2 + /«(x, 0) 2 >[po + pi//(l + r(r 1 //))] 



and in the region r > a, (i = 2) Laplace's equation 
holds exactly. This is the critical region for heat loss. 



(A-8a) 
(A-8b) 



292-251 0-68— 2 
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Now C is assumed small by virtue of 7 2 /(/f — J%) and 
z\(x) being small, so that 

BCg<$\Bg. 



Also since 



Pig^Po + PiG, 

BCg<AC. 

Thus we can neglect the second source term. We 
shall also neglect &(I 2 z\(x)) in C. Then we have 

k i0 ^g = -8n{g^ 2 Pi(l^r(r l G) 

+ r 1 GT'(riG)) + [po+piG(l + r(r 1 G))]- 

[/ 2 +(^->f)2i(*)]} (A" 11 ) 

We shaU deal only with G(A-8) and g(A-ll) to 
which must be appended boundary conditions: 

r= : 3g=3g = 0, (Q^x^l): 

dr dr 

l=b:G = g=0, (0«£*ssZ); 
x = 0.l: G = g=0, (O^r^b); (A-12) 

while the interface continuity conditions become 



tion at worst, we assume that the average axial value 
of G may be used to evaluate r(rjG): 



T(nG) = T (nG) (A-16) 

Then we shall use the differential equations 

ik i0 ^G + J^jo, ( 1 + t, ) 8,-jC = - J^?p 8a 

• [po + p.GU, 0)(1 + t,)]o h (A-17) 

where Ti = r(riG), t% = T\ + r\Gr' (nG) . 

On r= a, we have the interface condition: 



: a: 



dG 
'dr 



kJ-% 
dr 



0, 



= 0, 



(A-13) 



for continuity of heat flow, and 

r=a:|G(l + T(nG))|g^ = 0, 

\g(l + T(r&) + r&r'(r£))\£ 



0, 



(A-14) 



where we have again used \g\ < \G\ in the second 
condition in (A-14). 

Now even in the region of the wire where the tem- 
peratures are highest, r(riG) —0.1, so that the non- 
linearity is small. Furthermore, because of the high 
conductivity of the wire compared to the air, the 
radial temperature gradient in the wire is negligible, 
and 



G(r, x) = G(0, x) 



(A-15) 



inside the wire. 

The equations are still nonlinear in a way prevent- 
ing the easy use of linear methods of solution. Extend- 
ing the assumption that r(riG) is only a small correc- 



\G(a~, x) = G{a + , x)cr\ 



^here 



(Ti = 



a 2 = - 



\g(a-, x)=g(a + , x)ar 2 , 
l + r(r 2 G*) 

1 + Ti 

T- r(r 2 G*) + r 2 GV(;fG*) 

1+T 2 



(A-18) 



where G* is the number which gives continuity of 
temperature at r— a + in the mean by (A— 7a): 

G*(H-r(r 2 G*))=G(l + r(r 1 G)). (A-19) 

The system of equations we solve then, (A-17), 
subject^ to conditions (A-18, A-13, A-12), are linear 
when G is given. From the solution, a new G is com- 
puted, and the solution iterated to make G stationary. 
With a reasonable_guess, say only 25 percent off, five 
iterations brought G to within 1 percent. 

No attempt has been made to bound the uncertainty 
of the substitution error due to the assumptions made. 

Possibly a better choice for the constant G than that 
given in (A-16) might be made. If a 5 percent increase 
in G is used to compute the (r/, a,), the substitution 
error changes by 0.06 percent in a typical case. The 
approximation is clearly better for large l/b, for which 
the temperature profile is flat for most of the wire 
length. 

The constant / in the source of (A-17) for g is deter- 
mined by the requirement that the total bolometer 
element resistance be unchanged when the d-c current 
J'i is replaced by the current ((J'2^~^'k(x, t))). 
Equating these resistances: 

r^[po + PiG(0, *0(l+T<riG(0, *)))] 

Jo 

= f 1 dx[p + pi//(0, x)(l + T(rx//(0, *)))] (A-20) 
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Expanding //to 1st order in g, and again assuming 
G(0, x) is flat enough to be replaeed by its average in 
t(tiG), we find 



C0>(*) = -(1 + Tj )p,.^ f' ^ f"d V r,x ( Hx, 0; f, ij). 

Jo Jo 

(A-28) 



/: 



<fog(0, *)=0; 



(A-21) 



We define the numbers 



for the balanced bridge this relation determines 
/ in (A-17). 

When no rf power is introduced, the d-c power 
developed in the bolometer element is 



•■-/' 

Jo 



dx(\ + &*\x)) 



Wn 



& 2 >(x)_ 
1,2= T dx(l + G^\x))z x (x) 

= *naW\ dxlpo + pitGiO, x)(l + t^G))] (A-22) The bridge balance condition (A-21) leads to 

Jo 



] 

G* 2) (*) 



A-29) 



When the d-c current is reduced to J 2 and rf current 
is added to rebalance the bridge, the total power 
developed is 



P 



B 2 



J\-Jl A 2 



(A-30) 



77Y/ 2 ^[/ 2 2 -fi^ () 2 2 W] 



W T 



• [po + pi//(0, *)(1 + T(rifl))] (A-23) 

-r DC2 + r rf 

where by virtue of (A— 20), 

r DC 2=s 2 r DC „ t-y^i. (A-24) 

The substituted power is defined as 

J^sub = ^i)c - ^dc2 = (1 - C 2 Ww (A-25) 

The change in power developed AJT (with A, B, C 
as before (A-10), and using the same approximations 
and the bridge balance relation (A-21)) is found to be: 

aw = w T0T - Wnci = w n - r SUB 

= ira 2 I' dx{AC + J?\Bg) 
Jo 

= 7ra 2 | cb[jS) + piC(O f x)(l + Ti)] 

■[P+^f-^Izi^)]- (A-26) 
Then the substitution error is 



In terms of the (Au #/), the power terms can be writ- 
ten as 



AJT=7ra 2 p yf(l-p 2 ) 



^sub =(l-p 2 )7ra2p ^ 2 ^ 

/ 2 



j^ 2 - 4 2 



/4i+fli 



so that the substitution error is 



A l B>-B l A 2 
AtAo 



AtBz-BtA* 



AtAo 



(A-31) 
(A-32) 

l - (A-33) 



The axial temperature distributions themselves are 
found from 



IITTX 



~ J n=l 



AT(x) = T n (x)-T x (x) 
_ 7rp Q 



A17TX 



m^^T-l 1 -! 



.£' 



(A-34a) 



(A-34b) 



£■ = - 



AT 






i+ 



A^ 



(A-27) 



where £* is a sum over odd rc only, the coefficients 
<p$ are given in Appendix B, and 

= ~f- (B 3n A 2 -A 3n B 2 ), (A-35) 



Our primary concern is with a sinusoidal current 
distribution 



Let x (j) (*, r ; £, V) De Green's functions for the sys- 
tems (A-17) with the superscript (/= 1, 2) referring to 
(ti, o-i) for G and (r 2 , o- 2 ) for g. The Green's function is 
symmetric in (x, f ). Let: 



z(x)= cos 



n H>- 



(A-37) 



from which we compute z\{x) by (A-9a). 
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Then 



z 2 (x) = cos 2 <p cos 2 kl (- — -t 

+ sin 2 <p sin 2 jfcZ (|-f) + ^M, (A-38) 



^here 



J^(X): 



1 . 

2 



sin 2<p sin 2A:/ (~ — -7) 



does not contribute since it is anti-symmetric on (0, /) 

while G^(x) are symmetric. 

We find, averaging and noting &tf(x)=0. 

z 2 (x)= cos 2 cp 



. kl kl' 

1 Sin 2~ C ° S Y 
2 + A/ 



H- sin 2 <p 



. A/ if 
sin — cos — 

2~ 



A/ 



(A-39) 



We obtain Zl (x) from (A-38, A-39) using (A- 9b). 

In Appendix B, details of the calculation of the 
G ii] (x) by Fourier series, description of a program used 
to compute the (Ai, Bj) and the temperature distribu- 
tions are presented. 

Appendix B. Green's Function 

We are to determine the Green's function 
X {j \x, H £, 77) satisfying the system 

kto&x^ + S\ Pi (1 + Tj) X ^ = - S(x - f) • 8(r - r)) • 8a 

r (B-l) 

subject to the conditions 
X (J) (x, r; £, tj)=0, 
dx U) 



r=b, OsSxsS/; 
x = (0, 1), O^r^b, 



dr 



(*,0;{,ij)=0, 



(B-2) 
(B-3) 



,0 M 



,r-f-i -fa.W 

\ dr / r=a - \ dr / r=a + 



which satisfies the conditions (B-2) at 1= (0, /). We 
obtain the equations (dropping the superscript (/)): 

*w " 4 r ^ + (J^?pi(l + t) - £ioa 2 ) X r 



= sin a£ 

k i^ r ^ T 

#20 1 r 1 
r ar ar 



8(r— »?) 
1? 



(r < a); 
feoa 2 x r = 0, (r>a). 



(B-7a) 
(B-7b) 



For r>a, the solution satisfying (A-2) is: 
X T (a, r) = D, 



/o(ar)-^-K (ar) 
A.o(a6) 



(B-8) 



For s£ r < 17, the solution satisfying (A-3) is: 

X T (a, r) = D Jo (ye), (B-9) 

where 

e^^-Ca/) 2 , 
/B*«.*5;pifl(l+T)/*io. 

For 17 < r < a, the general solution is: 
x T (a, r) = DJ (je\ + DJ (j< 



(B-10) 



(B-ll) 



The four constants Z); are readily found from con- 
ditions (B-4, B-5) at r = a, and from jump conditions 
at r = n: 



X T (a, V )=X T ( a , V + )> 

d X r ( _._ a X r / + x sin qg 
— - (a, 17 )- dr («» ^ ) ZI 



We obtain for the axial value 



(B-12) 



(X (j) )r=a-=Cr,(xO ) )r=a + . 



(B-4) x(0 , r; £, ,) =7 f- J 

(B-5) 



R« 



A Fourier sine series for the solution is obtained by 
taking the finite Fourier transform. 

X T (a, r; £, 17)= c?xx(x, r; £, tj) sin ax, 
Jo 

2 °° 
X(x, r; f, i7) = y ]T x r (a«, r; f, 17), 



yo(^6„) + y () (^e 



. 7TTl£ . 7rnx 
sin — -r 3 - sin — — , 



(B-13) 



where 



R,= 



(f«.)c; + ^j, 



• (? e ") 



*(f*)« + S&*(7*) 



T" 



(B-6) 



„g _ Io(aan)Ko(ba„) — Io(ba H )Ko(aan) 
" /i(aa„)Ko(6a») — 7o(6a H )K r (aa„) 



(B-14) 
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All that is required in the text is the integral of this 
on 7). We obtain 



f° <W%> 0; ft t») = ^4 £ "*P -" =J* sin SJ 

Jo ^KIOPJ ~ « l 



(B-15) 



*g>- 



7Ttt€t 



'•(? 



tin + 



&j€jnki 

nirk'zo 



r,(f«.) _ 



(B-16) 



Note that the integral (B-15) is symmetric in (x, £). 
From (A-15), we obtain 

G»(*) --pi/Kl + ij) f d| f dnxf»{x, 0; ft ij) 

Jo J<> 

= -Pi^(1 + tj) [ <# J Awx^tf, 0; *, 17) 

= jj*$/>sin^ (B-17) 



The Bessel functions occurring in 0j^ are computed 
from polynomial approximations [9]. 
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